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A photon breeding mechanism for the high-energy emission of 
relativistic jets 
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ABSTRACT 

We propose a straightforward and efficient mechanism for the high-energy emission of rel- 
ativistic astrophysical jets associated with an exchange of interacting high-energy photons 
between the jet and the external environment. Physical processes playing the main role in 
this mechanism are electron-positron pair production by photons and the inverse Compton 
scattering. This scenario has been studied analytically as well as with numerical simulations 
demonstrating that a relativistic jet (with the Lorentz factor larger than 3^1) moving through 
the sufficiently dense, soft radiation field inevitably undergoes transformation into a luminous 
state. The process has a supercritical character: the high-energy photons breed exponentially 
being fed directly by the bulk kinetic energy of the jet. Eventually particles feed back on the 
fluid dynamics and the jet partially decelerates. As a result, a significant fraction (at least 20 
per cent) of the jet kinetic energy is converted into radiation mainly in the Me V - GeV energy 
range. The mechanism maybe responsible for the bulk of the emission of relativistic jets in 
active galactic nuclei, microquasars and gamma-ray bursts. 

Key words: acceleration of particles - galaxies: active - gamma-rays: bursts - instabilities - 
methods: numerical - radiation mechanisms: nonthermal - shock waves 



1 INTRODUCTION 

Traditionally, the dissipation of a relativistic bulk motion into 
radiation is assumed to be associated with the shock accelera- 
tion of charged particles. In a relativist ic case its role is lim- 
ited by a numb er o f facto rs (see e.g. JA chterberg et al. 2001; 
Bednarz & Ostrowskil Il999l) . iDerishev et alj J2003l) "and ISternl 
(2003) independently suggested that interacting neutral particles 
can convert bulk kinetic energy into radiation much more efficiently 
than this can be done by charge particles. Indeed, the neutral par- 
ticles easily cross the shock front or the boundary of the shear 
layer in both directions and can be converted into charge particles 
for example via e pair productions by two photons. Pairs in turn 
convert their energy into photons by Compton scattering, some of 
which cross the boundary again. Such exchange of particles be- 
tween media moving respectively to each other with high Lorentz 
factor works in the same way as Fermi acceleration, but more ef- 
ficiently as there is no problem with diffusion through the bound- 
ary layer. This process can proceed in a runaway manner where 
the high-energy photons breed exponentially similarly to neutron 
breeding in nuclear chain reaction. 

Stern 1 2003) demonstrated with numerical simulations that in 
the case of an ultra-relativistic shock in a moderately dense medium 
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(above ~ 10 4 particles cm -3 ) this process leads to a dramatic in- 
crease in the shock high-energy emission (electromagnetic catas- 
trophe) and eventually to the elimination of the shock front, which 
is converted into a smooth radiation front. Such effect can take 
place in gamma-ray bursts. In this work we consider the same 
scenario for a different kind of relativistic fluid: a shear flow in 
astrophysical jets. While we specifically consider here jets in ac- 
tive galactic nuclei (AGNs) which emit tremendous power in form 
of the gamma-rays in the MeV-TeV energy range, the mechanism 
can operate in the jets from microquasars as well as in gamma-ray 
bursts. 



There exist different models of the jet emission mechanism, 
probably th e most popular one is associated w i th internal s hocks 
in the jet lReedll978t IPaczvnski & Xulll994 iRees & Meszarod 
1994). Nevertheless, the maximal variation of V is at the jet bound- 
ary and it is there we can expect the most intensive energy re- 
lease. The jet bulk energy can be dissipat ed in principle th rough 
the a c celeration of charged parti cles (Berezhko 1990c lOstrowskil 
2000; Stawarz & Ostrowski 2002). If Thomson opacity is suffi- 
ciently large, the standard radiation viscosity can also cause the 
energy dissipatio n in the bou ndary layer as, for example, discussed 
bv lArav & BegelmanH l992). In the case of AGN jets the Thomson 
opacity is, however, insufficient to provide the efficient reflection of 
photons crossing the boundary of the jet. Here we show that the role 
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of the "mirror" can be played by the soft radiation, which provides 
the pair production opacity for the high-energy photons. 

We study the mechanism of shear flow energy dissipation nu- 
merically in the same way as it was done by Stern 1 2003) for 
shocks. The problem formulation differs in its geometry (tangen- 
tial jump of the fluid velocity instead of the head-on shock) and in 
the character of the primary photon background. The complete so- 
lution of the problem should account for the feedback of particles 
on the fluid dynamics and requires detailed numerical treatment 
of the hydrodynamical part of the problem. In this first study we 
neglect internal pressure and consider a one-dimensional dust ap- 
proximation. With this simplified treatment of hydrodynamics we 
formulate our problem in the following way. We start from a non- 
radiating, idealized cylindrical, homogeneous jet and then track its 
evolution to check whether it produces a runaway photon breeding 
or not. This numerical experiment is sufficient to demonstrate that 
under certain circumstances a non-radiating jet cannot exist: it must 
radiate away a substantial fraction of its kinetic energy. 

We present the model of the jet and the environment and for- 
mulate the quantitative problem to be solved in Section|2| In Sec- 
tion[3] we describe the proposed mechanism qualitatively, formu- 
late the necessary conditions needed for its operation, and provide 
a simplified analytical study of separate parts of the mechanism. 
In Section|4| we give the details of the numerical method for sim- 
ulation of the entire process. The results of numerical simulations 
for two representative sets of parameters, which show a runaway 
regime are presented in Section|5| We discuss the results, problems 
with the model and possible effects in more realistic situations in 
Section|6| and we conclude in Section^ 



2 MODEL OF THE JET AND THE ENVIRONMENT 

We consider a jet of Lorentz factor T consisting of cold electrons 
and protons and the magnetic field. We consider a piece of the jet 
centered at distance R from the central source and length 20Rj, 
where the jet radius is Rj = R6 and 8 has the meaning of the 
opening angle, while for simplicity we approximate the jet by a 
cylinder of radius Rj . We choose parameters corresponding to the 
case of AGN taking the distance scale R ~ 10 17 cm, T ~ 10- 
20 and = 0.05. The radial distance from the jet axis r and the 
distance from the central source z is measured in units of Rj . The 
unit of time is Rj/c. 

The jet propagates through the soft radiation field, which en- 
ergy flux we denote F(x), so that F(x)dx is the energy flux in the 
interval da;. The photon energies, in units of the electron rest mass 
m c c 2 , are denoted as x and e for photons of low (< 1) and high 
( > 1) energies, respectively. We define also the power-law index 
a = — dlogi^a^/dlogi. 

The radiation field in AGNs has two major components. The 
first component originates in the accretion disc around the cen- 
tral source and its photons propagate along the jet direction. We 
consider a standard accretion disc iShakura & Sunvaev 1973), as- 
suming a simple power-law T(R) ex R~ 3 ' 4 dependence of tem- 
perature on radius, with the ratio of the outer to inner disc radius 
-Rout /Rin = 10 4 and the maximum temperature of T = 5 eV. The 
dimensionless temperature is then = kT/m c c? fn 10~ 5 . The 
resulting multicolor disc spectrum Fd(x) has a power-law shape 
with energy index —1/3 below B and has a Rayleigh- Jeans part at 
energies below e min = 6 (Rout/Rin)~ 3/4 ~ 10" 8 . 

The second, isotropic component consists of the disc radiation 
scattered and reprocessed in the broad line region (BLR) and the 




Figure 1. The scheme of the photon breeding cycle. The cylindrical jet 
moves to the right with Lorentz factor I\ The size of the Larmor orbit is 
highly exaggerated. See Sectionl3lfor details. 



infrared radiation by the dust JSikoraet alJll994 Il996h . We take 
its energy density to be a fraction r\\ of the energy density of the 
direct disc radiation, and its flux to be a cutoff power-law F\ so (x) oc 
x~ a exp(— i/i min ) extending from the far-infrared aw n ~ 10~ 9 
to the UV band, as we assume a; max = B. 

The jet kinetic luminosity scales with the disc luminosity as 

2 



<r) K L d = T Mc* 



(1) 



From the mass conservation (for a two-sided jet) M — 
2tt8 2 R 2 rn p cn p (R), we get the Thomson optical depth across the 
jet from the electrons associated with protons 



ttOR) = n p (R)a T Rj = 2.4 KrVL d ,45/(fli70r), 



(2) 



where we use standard notations Q — lO^Q^ in cgs units. 

Comoving value of the magnetic field in the jet is Bj, its di- 
rection is transversal by assumption. A reasonable scaling for the 
jet field is Bj oc 1/Rj, then the Poynting flux carried by the jet 
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x W 43 B 2 R 2 17 (eT) 2 erg s _1 (3) 



Lb = mLd = -±2wRJT 2 c ■ 

on 

is constant along the jet. Parameter r/B defines the role of the mag- 
netic field and we consider two cases with low and high tjb ■ 

For our analysis we consider two reference frames: the am- 
bient medium (external) frame and the frame comoving with the 
fluid with Lorentz factor T (jet frame). Energies are normally given 
in the external frame. The jet comoving values are specified using 
subscript 'c' or a prime. 



3 STEP-BY-STEP DESCRIPTION OF THE BREEDING 
CYCLE 

The solution of the problem of the jet interaction with the external 
environment and dissipation of its kinetic energy into high-energy 
radiation requires full-scale nonlinear Monte-Carlo simulations as 
described in Section [4] However, it would be useful to perform a 
simplified analytical study. Very schematically, the process can be 
split into five steps (see Fig.0. 

(1) A high-energy external photon (which origin is not impor- 
tant) enters the jet and interacts with a soft photon producing an 
electron-positron pair. 

(2) In the jet frame, the produced e + and e~ originally move 
backwards relative to the direction of jet propagation with the en- 
ergy ~ F times higher than the energy of the parent photon (in the 
external frame). When particles gyrate in the magnetic field of the 



Photon breeding in relativistic astrophysical jets 3 



jet, their time-averaged energy, as measured in the external frame, 
becomes ~ Y 2 higher than the parent photon energy. 

(3) The pair Comptonizes soft photons up to high energies. 

(4) Some of these photons leave the jet and produce pairs in 
the external environment. 

(5) Pairs gyrate in the external environment and Comptonize 
soft photons more or less isotropically. Some of the Comptonized 
high-energy photons enter the jet again. This step completes the 
cycle. 

Steps 1, 3, 4 and 5 require a field of soft photons to provide 
the conversion of the high-energy photons into pairs and to produce 
new high-energy photons through the inverse Compton scattering. 
In our case, the direct disc radiation provides the opacity for the 
high-energy photons at step 1 and provides a target for Compton 
scattering at step 5. Because the disc radiation is directed along the 
jet, its effect at steps 3 and 4 is small (interactions are suppressed 
by the factor ~ 1/T 2 ). The isotropic soft photon component pro- 
vides the major source of seed soft photons for Comptonization at 
step 3, provides opacity at step 4 as well as contributes to the pair 
production at step 1 . 

An additional condition for high efficiency of the cycle is the 
presence of a transversal or chaotic magnetic field, both in the jet 
(to provide step 2) and in the external environment (to provide 
isotropization of pairs and photons Comptonized by them at step 5). 

Each step can be characterized by its energy transmission co- 
efficient Ci defined (in the external frame) as the average ratio of to- 
tal energy of particles with energy above the pair production thresh- 
old before and after the step, c-2 is large (~ T 2 ), others are smaller 
than 1. If the criticality index C = ci x C2 x C3 x C4 x C5 > 1, then 
the regime is supercritical, i.e. each cycle produces more particles 
than the previous one and their number grows exponentially. In this 
case we deal with particle breeding rather with particle accelera- 
tion. The spectrum of particles changes slowly (and in principle the 
mean particle energy can decrease), but the number and the total 
energy of particles grow rapidly. 

Each coefficient c k can be represented as an average of energy 
dependent coefficient C k (e) over the photon spectrum at fc-th step: 



Ck 



C k (e)F k (e)de, 



(4) 



where F k (e) is the energy flux spectrum at fc-th step normalized 
to unity. The spectrum is a matter of a self-consistent treatment of 
the whole breeding cycle and we cannot estimate it a priori. For 
simplicity, we estimate C*.(e) for the constant external photon field 
(disc as well as isotropic). 

A quantitative example for Cfe(e) behaviour described in this 
section is given for the same parameters as in the detailed simu- 
lation presented in Section BTl id, 45 = 1, R17 — 2, T — 10, 
= 0.05, tjk — 1, TIb — 0.01, 771 = 0.05, and we consider here 
two cases of the isotropic photon spectrum with a — 0.4 and 1. 
The value a = 0.4 seems reasonable for the BLR region, where 
the scattered radiation should be dominated by the UV and optical 
component and some less energetic IR radiation of dust is expected, 
while a = 1 could be reasonable at the parsec scale where dust ra- 
diation can be very important. 



3.1 Steps 1 and 4: propagation of photons and pair 
production 

Let us consider a problem of photon propagation through the 
boundary between the external medium and the jet. Let a 77 (e) be 
the total opacity for pair production (i.e. the absorption coefficient 



produced by the disc and isotropic components) for a photon of 
energy e. The opacity produced by the isotropic component of the 
radiation field only is denoted a 77 (e). 

The jet photons move within the cone of opening angle 1/T to 
the jet direction and interact mostly with the isotropic component 
of the external radiation field, because the probability of interaction 
with the disc component is reduced by the factor ~ 1/T 2 . The cor- 
responding opacity computed for the interval perpendicular to the 
jet direction is <jj ~ ra 77 (ej). The probability that the jet photon 
is absorbed in the interval dzu at distance vj from the jet bound- 
ary is exponential exp(— ajza)a,jdzci (this assumes that all photons 
cross the boundary at the same angle 1/T). The pair production and 
consequent Compton scattering transform this photon to photons 
of energy e c . Roughly half of these photons are emitted towards 
the jet and the fraction exp(— a c vj) reaches the boundary (where 
a c = a 77 (e e ) is the absorption opacity and we assumed that the 
photons propagate in the direction perpendicular to the boundary). 
The probability that the photons are absorbed within the jet is about 
1 — exp(— a c Rj). Assuming that only photons born at a distance 
not more than the jet radius from the boundary can enter the jet, we 
get the total probability for an external photon of energy e c to inter- 
act inside the jet with a soft photon producing an electron-positron 
pair: 



e J a\dzu 1 — e ' 



<"L(e e ,ej) ~ / e °- ro | e a J ro aj dw [ 

- iUsi 1 --*]! 1 - 



-(a c + aj)fij 



where 

Oj_ _ ra 77 (ej) 



(e.) 



, (5) 



(6) 



Obviously, C\ depends also on the energy of the parent photon q , 
because it defines the spacial distribution of sources of photons e c . 
Analogously we derive the probability of a high-energy pho- 
ton ej produced in the jet to escape from the jet and to produce 
a pair in the external medium at a distance not more than the jet 
radius from the boundary: 



C 4 (ej,£c) 



z ' a c dvj 1 



lihV-^W 



[l_ e - a J^] 



-(a e +aj)iij 



(7) 



Again, C4 depends on the energy of the parent photon e c that gave 
rise to the photon e- s in the jet. 

If opacity is sufficiently high (i.e. we can ignore exponen- 
tial factors in equations J5J and 0), the sum of the probabilities 
Ci(e e , £j) + C4(ej, e e ) = 1/2, and their product reaches the max- 
imum of 1/16 when the opacities are equal, b — 1. Accounting for 
the angular distribution of photons gives 



max{Ci(ee,£j) C4(ej,e e )} « 1/21. 



(8) 



Let us now discuss in details the behaviour of the photon opac- 
ity. For isotropic external photons, the main contribution to the pho- 
ton opacity is given by the radiation of the accretion disc F<x(x). 
The angle-averaged cross-section of photon-photon pair produc- 
tion cr 7 -y(e, x) has a low-energy threshold x e > 1 and reaches the 
maximum « 0.21<jt at xe sa 3.5 (see e.g. Zdziarski 1988). The 
absorption coefficient for a photon of energy e can be obtained by 
integrating the cross-section over the spectrum of soft photons: 



a 77 (e) = A r phO- T s 77 (e), 



(9) 
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Figure 2. The averaged over the photon distribution photon-photon absorp- 
tion cross-section s 77 (e) given by equation 1101 . The dotted curve gives the 
cross-section averaged over the multicolor disc spectrum with © = 10~ 5 
and R ut/R\n = W 4 . Almost straight lines give the opacity induced by 
the isotropic component of the radiation field oc e a exp(— l/e0) (scaled 
by the ratio of photon densities A r p i 1 4 so /-'Vph,d of the isotropic and disc 
radiation). Thick curves give the sum of the two opacities. Solid and dashed 
lines are for a = 0.4 and a = 1, respectively. 



where the mean cross-section (in units of or) 
s 77 (e) = / n(x)a 1J (xe)dx 

Jl/e 



(10) 



and n(x) is the photon number density normalized to unity. The 
mean cross-section for the multicolor disc is plotted in Fig. [3] by 
the dotted curve. The total disc photon number density can be 
expressed through the energy density iV p h,d = (JJa/kT, where 
£ ~ 1.4 for the multicolor disc. Because Ud — 2La/(c4irR 2 ) 
close to the jet axis (factor 2 comes from the angular distribution 
of the disc radiation, which we assume to follow the Lambert law), 
we obtain 



,(0 



(TT 1 

cA-kB? m e c 2 Q 

id, 45 



6 10" 



2£s 77 (e) 
-s 77 (e) cm~ . 



(11) 



The opacity is low at e < 1/0 ~ 10 and has a maximum at 
£ma X ~ 10/9: a 77 , max = 6 10~ 15 Ld,45/(-R? 7 e-5) cm -1 . At 
higher energies e > e ma x, the photons interact with the power-law 
part of the multicolor disc spectrum resulting in a power-law decay 
of the opacity: 



r(c) 



1/3 



(12) 



which transforms into a faster - In e decay at e > l/O m i n - 

In spite of a lower energy density, the isotropic component 
Fiao(x) dominates the opacity at very high e, because of its much 
softer spectrum a > —1/3. The opacity increases with energy fol- 
lowing a power-law (Svensson 1987) 



W - i 



1 0"T 



5 m e c? 



Ftao(l/c) 



(13) 



as shown in Fig. [2] 

Because the coefficients C\ and d depend on both pho- 
ton energies e c and q, it is easier to visualize them as one- 
dimensional cross-sections. As an example, we plot Ci(e e , 10e o ) 
and C4(ej,ej/10) in Fig. [5] Coefficient C\ vanishes below e c <C 
E m „, because of low opacity when the photons escape freely and 
do not produces pairs within the jet. In the interval 5 10 4 < e c < 
10 6 , C\ ~ &/2(l + b) ~ 1/4 because here the ratio of opaci- 
ties b ~ 1. At higher energies, the opacity is dominated by the 
isotropic component, b ~ ra 77 (ej)/a 7 °(e G ) ~ r(ej/e c ) a is large 
and Ci ~ 1/2. 

Coefficient C4 vanishes at q < 2 10 5 because the opacity 
produced by the isotropic power-law component is low here, and 
photons are escaping too far from the jet (note the exponential fac- 
tor in eq. Q). At q ~ 10 6 , b ~ 1 and C 4 ~ 1/4. At high en- 
ergies d decays oc l/2b. In our formulation of C4, we accounted 
only for photons that can directly penetrate from the jet to the exter- 
nal medium. However, those high-energy photons that are absorbed 
within the jet do not leave the system, but produce pair-photon cas- 
cade. During the cascade the photon energy eventually becomes 
sufficiently low to allow the photon escape from the jet and to pro- 
duce a pair outside. Thus even if a photon was emitted towards the 
jet axis, a fraction of its energy can eventually escape the jet in a 
form of secondary photons. 

3.2 Step 2: energy gain of a pair produced in the jet 

A high-energy photon of energy e c at this stage has interacted with 
a soft photon inside the jet to produce an electron-positron pair. The 
two particles move now backwards relative to the jet propagation 
direction with the mean comoving Lorentz factor 



7< 



= — -r. 
3 2 



(14) 



Here the factor 4/3 comes from the averaging over the angles of 
the incident high-energy photons accounting for the probability of 
interaction with the disc photons moving along the jet (in case of 
interaction with the isotropic photon field, this factor disappears). 
Integrating over the Larmor orbit, we get the mean particle energy 
in the external frame 



< 7 > = 7 J {(1 + cos0 c ) 2 ) 0c = ^ 7c r = 2|r 2 , 



(15) 



where <j> c is the angle between the particle momentum and jet prop- 
agation direction (in jet frame). Note, that the averaging is over c/> c , 
not cos 4> c ■ The mean energy gain is then 



C 2 (e B 



_(7>_ 
e c /2 



= 21^ 



(16) 



The time-scale of synchrotron cooling is equal to the inverse 
of Larmor frequency at the (comoving) electron Lorentz factor 



7c 



10 S R 



-1/2 



(17) 



Compton losses because of the disc radiation can be neglected at 
this step (for our parameters) due to a deep Klein-Nishina regime 
of scattering at j c > 10 8 . For very high initial photon energy e c 2> 
7c,ma x r, the Lorentz factor of the produced pair is reduced to ~ 
7c, max after first Larmor orbit because of the synchrotron losses. 
In the external frame its energy is now ~ 7c,maxr, resulting in the 
energy gain (loss) factor C2(ec) « 7c,ma x r/e c . Thus the energy 
gain for arbitrary e c (shown in Fig.|5J can be written as 

C 2 (e c ) = 2r 2 T — — i;- . (18) 



1 



=r/ 7 c 
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where 



Figure 3. Coefficients Cj(e) versus photon energy. C2 and C3 are func- 
tions of energy of a photon produced in the external medium e e , while C5 
is the function of a typically higher photon energy e- y produced in the jet. 
Coefficients C\ and C4 depend on both energies, and we plot here only 
the one-dimensional cross-sections C'i(e c , 10e c ) and Ci{e- lt ej/10). Solid 
curves represents the case a = 0.4, while dashed curves are for a = 1 (C2 
and C5 do not depend on the choice of a). 

3.3 Step 3: energy conversion into photons 

Coefficient C?,{e c ) describes a fraction of the electron (positron) 
energy, which is converted into photons with the energy above the 
thresholds of the cycle £j, m i n = 3 10 5 defined by a sharp drop 
of coefficient C4. Photons of lower energies escape almost freely 
from the jet and are lost from the cycle. C4 also drops above ~ 
3 10 7 , however, the photons of higher energies are not lost from 
the cycle, but produce a pair cascade and eventually escape from 
the jet producing pairs outside. 

The efficiency of conversion is mainly determined by the ra- 
tio between Compton and synchrotron losses (because synchrotron 
photon energy is below £j, m m in our case), which, in Thomson 
regime, is just the ratio of the soft radiation energy density (dom- 
inated by the isotropic component) to the magnetic field energy 
density (measured in the jet frame): 



mu d r 2 



UL 



U^ L b /(2ttR 2 T 2 c 



= HL(re) 2 r 2 



VB 



(19) 



In reality, however, the share of Compton scattering is much 
smaller, because a pair produced in the jet by a photon at e c ~ I/O 
gains energy up to 7 ~ T 2 /S and interacts with photons at x ~ Q 
in the deep Klein-Nishina regime. The situation also complicates 
by the fact, that because of cooling, the fraction of soft photons 
interacting with the pair in Thomson regime grows with time. 

Neglecting for simplicity scattering in the Klein-Nishina 
regime and pair cooling (i.e. assuming constant pair Lorentz fac- 
tor given by Eq. 1141 ). we define coefficient C3 as the ratio of the 
soft radiation energy density, that can produce photons above the 
threshold, to the sum of the magnetic field energy density and the 
total radiation density in the Thomson regime: 



C 3 (ee 



Tjl 

w cycle 

U'+U' 



(20) 



U B = B 2 /(8n), 

f^KN 

Ut cx T 2 / F iBO {x)dx, 



^cycle ^ ^ 



Fis (x)dx. 



(21) 

(22) 

(23) 



maxhr,..;,. ,xi. 



The upper integration limits is defined by the Thomson regime of 
scattering 



x < xkn 



1 



(7) ^ 



(24) 



The lower limit is given by the condition that the scattered photon 
energy q = ^(-y 2 ) is above the threshold ej, m i n . Averaging over 
the Larmor orbit, we get the mean square of the particle Lorentz 
factor: 

(7 2 > = 7?r 2 ((1 + cos C ) 3 )„ C = | 7c 2 r 2 = H e 2 r 4 , (2 5) 

that gives 



X J> Xi ow — 



9 q,. 



<7 2 ) 



10 e 2 T 4 



(26) 



The behaviour of Cs(e e ) is shown in Fig. [3] 

There are three characteristic energies. For e c < ei ~ 
ei.min/r 2 ~ 10 3 , ikn < x\ „ and C3 vanishes. At energies 

10 5 , 



£c > £2 



'low 

^2 



\/ e J,miii/2;min/r 2 ~ 10°, Ut = [/"cycle because 
xi ow < aWn. And at energies above £3 = l/(T 2 a; m i n ) ~ 10 7 , 
ikn < Kmin and C3 = in our approximation. At energies 
£2 < e c < £3, one can approximate 

/ v, r — r 

1 + 



1 



C 3 (£e 



•' -11 



1 + 



uk_ 

Tjl _l-a _ 1-a ' 

^rad X KN ^min 

^B m(^max/^min) 



Q#l, 



(27) 



C/ r ' ad ln(a;KN/a;min) ' 

Taking into account equation 1191 . we can conclude that 
Comptonization at step 3 can be very efficient (C3 ~ 1) if: 

(i) The magnetic field is weak, ?Jb <S 1. 

(ii) The energy density of the external isotropic radiation is large 
in the jet frame, i.e. V or rji are large. 

(iii) The spectrum of isotropic component is soft, a > 1. Then 
the power-law decay of C3 above £2 is replaced by a constant or a 
logarithmic decay (see eq. 1271 and notice a slow decays of C3 in 
Fig. [3] dashed curve). In this case, a high value of C3 is possible 
even for magnetically dominated jet if its total power is a few times 
smaller than the disc luminosity. 



3.4 Step 5: Compton scattering in the external environment 

At step (5) an electron-positron pair, produced in the external en- 
viroment, Comptonizes a number of soft photons to produce high- 
energy photons above the energy threshold £ e ,mm ~ 3 10 4 defined 
by the low-energy cutoff of C\. At this step, synchrotron losses 
can be safely neglected as magnetic field in the external medium is 
much lower than the jet field. 

Because cooling is fast, the spectrum produced by a pair is 



a standard "cooling" spectrum F(e) 



-1/2 



extending up to 



the maximum energy equal to the pair Lorentz factor 7 = Cj/2. 
Defining coefficient Cs(ej) as the energy fraction emitted above 
the thresholds £ c , m i n , we get 



Boris E. Stern and Juri Poutanen 



M 
o 



I I I I I I I l/l I I I I I I I I I I M I I I I I I I 




I I I I I I 



Figure 4. The product of the five coefficients A(e c ,€j) = 
Cl(e e , ej)C2(ee)C3(e e )C4(ej, e e )Cs(ej), represented by contours 
of constant levels. For a = 0.4 (solid contours), it reaches the maximum 
of 7.4 at e c = 2 10 5 , £j = 10 6 . In this case, the dotted curve shows 
the median energy £j(e c ) for steps 2 and 3, and the dashed curve shows 
the median energy ee(ej) at step 5. For a = 1 (dashed contours), the 
maximum of 9.9 is at e c = 3 10 5 , e; = 2 10 6 



n(> 



C S (€i) = 1 ~ 



1/2 



= 1- 



2fe 



1/2 



(28) 



The computed Cs(ej) is shown in Fig. [3] A low-energy cutoff is at 
ej = 2e c , min = 6 10 4 . 



3.5 Amplification through the cycle and dissipation efficiency 

The absolute theoretical maximum of the amplification factor 

maxyl(e c , ej ) = C1C2C3C4C5 « T 2 /10 (29) 

is achieved when C 2 = 2T 2 , Ci X C 4 = 1/21, and C 3 = C 5 = 1. 
Thus, the minimum jet Lorentz factor required to achieve super- 
criticality is T ~ 3. In a more realistic situation, the amplification 
factor is smaller and larger T is needed. For example, strong mag- 
netic field inhibits the photon breeding reducing C2 and C3 because 
of synchrotron losses. The cascade still develops, if a sufficiently 
dense soft external radiation field is present, because C3 ~ 1, when 
riir 4 /r)B is large (Section f3.3t , see Sect. l5.2l for the simulation ex- 
ample. 

An estimations of A using our simple formulae from Sec- 
tions |33031 i s shown in Fig. [4] For a = 0.4, it reaches 7.4 at 
e c = 2 10 5 , €j = 10 6 (the more precise Monte-Carlo step-by-step 
simulations of the cycle efficiency show A reaching the maximum 
of 4.7 at e e = 5 10 5 , ej = 2 10 6 ). For a = 1 the maximum of A 
is about 9.9 at e c = 3 10 5 , q = 2 10 6 . Of course, A(e c , 3) does 
not have the meaning of the criticality index C. The photon en- 
ergy distribution during the cycle is wide and does not necessarily 
coincide with the area of the maximal amplification. Fig. [4] (dot- 
ted curve) shows the median energy ej(e e ) computed by a Monte- 
Carlo method (for a = 0.4) of the distribution tjP23(ej, e e ), where 
P23(tj, e c ) is the probability density of a photon of energy e c to 
produce a photon of energy ej after steps 2 and 3. A similar func- 



tion e e (ej) for step 5, which is a median of distribution e c ps(e c , ej), 
is also shown (dashed curve). Note that the median energy gain at 
steps 2 and 3 is slightly higher than the energy loss at step 5, which 
means that photon energy rises on average during the cycle and 



tends to the area at e c 



3 10°, a 



3 10 , where the curves 



intersect. The amplification in that area is A(e e , ej) ~ 2. 

We can now estimate a typical distance where the most effi- 
cient energy dissipation in the shear flow should take place. The 
total optical depth across the jet is (see eq. 1111 ) 



r 77 (e) = a~ fl (e)Rj = 300 



L dA5 (2Q8) _ 
Ri7<d- 5 ' 



,(e) 



(30) 



which gives us the maximal distance from the central engine where 



r (tn 



1 and, therefore, is sufficient to produce pairs in the 



jet: 

L d , 45 (206) 

fe)-5 



(31) 



Note that infrared emission of the dust can provide the opacity for 
higher energy photons at much larger distances. 

At R < -Rmax the typical depth of photon penetration 5zo — 
8r Rj into the jet is given by the condition 5m a 11 (e) — 1, i.e. 
Sr = 1/t 77 (e). At the peak of absorption, for e ~ e ma x, we obtain 

R „„„ e_ 5 -Ri7 



5r = 



= 0.03- 



(32) 



Rmax Ld,45 (208) 

The value of 2Sr defines the fraction of the jet volume occupied by 
the "active layer", which is responsible for the exponential breeding 
of the high-energy photons. Because the photons of slightly higher 
and lower energies than e ma x can propagate further into the jet, the 
fraction of the total jet kinetic energy released is even higher that 
25r, as demonstrated by the full scale simulation (see Section |5}- 
At small distances from the central source, opacity is so large that 
high-energy photons can hardly penetrate the boundary between 
the jet and the external medium. The cascade can develop there, 
but the dissipation efficiency is very low because 8r is small. 

The efficiency also depends on the number of generations N 
the breeding cycle operates. A high efficiency is achieved if N is 
larger than a few (say N — 10). The time-scale of the cycle is de- 
fined by the longest step 4, which takes on average i cy cic ~ Sr T. 
Because the active layer width depends on distance Sr oc R, the 
breeding cycle takes longer time (in Rj/c units) at larger R and 
at the same time it is less sensitive to the requirement of a sharp 
boundary. The cycle time-scale should be compared to the dynam- 
ical time-scale tdyn ~ 1/0. This determines the most favourable 
distance for the breeding to operate as 

flmax 1 id, 45 20 



R 



Nve 



n e. 



pc. 



(33) 



For this estimation we assumed that the amplification factor 
A does not depend on distance R. This is true as long as B?, Fd 
and _Fi so all scale as 1/R? and Sr <g 1. The last condition breaks 
down at a parsec scale (see eq. 1311 ). while the scaling for F- lBO 
probably breaks much closer, as one can hardly expect the existence 
of the sufficient material to provide efficient scattering/reprocessing 
of the disc radiation to the isotropic component beyond the BLR. 
Thus, the most efficient dissipation is expected within the BLR at 

RcH < R < -Rmax- 

There still remains such source of _Fi so as the dust IR radiation, 
but it cannot convert photons at ej < 10 7 into pairs and the coeffi- 
cient d will be small. In such a case one can expect that the cycle 
works in the TeV range, which could be relevant for TeV blazars. 
At small distances R < 10 1B cm, there appears another possibility: 
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Figure 5. External frame (7, dashed curve) and jet comoving (7 C , solid line) 
Lorentz factor of a particle gyrating in transversal magnetic field Bj = 
0.35 G frozen into fluid. Upper panel: dependence on the external frame 
time. Lower panel: dependence on the comoving time multiplied by the bulk 
Lorentz factor T = 10. Initial Lorentz factor of the particle is -y c = 3 10 s 
and the jet radius Rj = 10 16 cm. 



nonthermal X-ray component of the disc radiation (which typically 
constitutes ~ 0.1 of the disc luminosity) becomes opaque and can 
efficiently convert lower energy photons (e < 1/6). Then one can 
expect the existence of a low-energy breeding cycle in the range 
e ~ 10 -10 . This possibility requires a separate study. 



4 NUMERICAL IMPLEMENTATION 

We assume a sharp boundary between the jet and the surrounding 
medium at the start of simulations. For simplicity we adopt the con- 
stant physical conditions along the 20i?j interval, including a uni- 
form field of the external radiation. The jet decelerates transferring 
its momentum to the radiation. In our model we take into account 
the dependence of the fluid Lorentz factor on r only. 

The numerical simulation method is based on the Large Par- 
ticle Monte-Carlo code (LPMC) developed by Stern ( 1985) and 
Stern 1 1995). The code can treat essentially nonlinear problems 
when the simulated particles constitute at the same time a target 
medium for other particles. The number of large particles (LPs) 
representing photons and e pairs was 2 17 = 131072. 

The version of LPMC used here treats Compton scattering, 
synchrotron radiation, photon-photon pair production, and pair an- 
nihilation. Synchrotron self-absorption was neglected as it con- 
sumes too much computing power and is not very important in this 
application. All these processes are reproduced without any simpli- 
fications at the micro-physics lev el. The gener al organization of the 
LPMC simulation is described in Stern] (1995J). A new specific fea- 
ture of numerical simulation in this work (as well as in Stern 2003) 
is a scheme of particle tracking in the relativistic fluid. 

Parameters of photon LPs are defined in the external reference 



frame. The energy and the direction of electron/positron LPs are 
defined in the jet frame, since the electron energy in the external 
frame oscillates by a factor 4L 2 because of gyration in the magnetic 
field. 

The tracking scheme for high-energy charged particle LPs dif- 
fers for the first Larmor orbit and the rest of the trajectory. The rea- 
son is that at the first orbit a particle can lose the main fraction of 
its (jet frame) energy before it turns around and gains the energy in 
the external frame (see Sect. 13. 21 and Eq. ll7i . 

Rapid energy losses require a fine particle tracking at the first 
Larmor orbit if the particle has been produced with 7 C > 7 c ,max- 
Therefore, the comoving tracking step is limited by As — 0.17?l, 
where Rl = 1.7 10 3 7 c /-Bj cm is the Larmor radius. Each step 
is described in both reference frames. The corresponding Lorentz 
transformations from the jet frame to the external frame are: 

dt = d£ c (l + — /3c cos 0c J T, 7 = 7c (l + — /?c cos 4> c J T,(34) 

where V is the velocity of the fluid, dt c is comoving time inter- 
val for the tracking step, /3 C is the particle comoving velocity in 
units of c, <f> c is the angle between comoving direction of the parti- 
cle momentum and the jet, At and 7 are the external frame values 
for the step time-interval and the particle Lorentz factor. The co- 
moving representation is used to track the particle gyration in the 
magnetic field frozen into the jet and to simulate the synchrotron 
radiation. The external frame representation is more convenient to 
simulate interactions with photons and is necessary to synchronize 
the particle tracking with the general evolution of the system. 

An example of high-energy {-y c > 7 c ,max) particle tracking 
for several Larmor orbits as viewed from both reference frames 
is shown in Fig. [5] The magnetic field and the external soft photon 
field correspond to the case considered in Sect. 15.11 One can see the 
dramatic energy loss due to synchrotron radiation at the first orbit 
and a much slower further evolution. A discrete step at t ~ 4 10 -3 
is a result of Compton scattering. 

When the particle energy is below 7 c , ma x, we can neglect the 
dependence between the comoving direction and the energy. In this 
case we sample the direction of the particle assuming a uniform 
distribution of its gyration phase <f) c in the comoving system. The 
external frame probability density function for </> c is 



p(4>c) oc dt/dtc 



1 + — /3c cos 4> c I r 



(35) 



Trajectories and momenta of LPs are three-dimensional. The 
target LP density is averaged over 45 two-dimensional cylindric 
cells: 5 layers along the jet with 9 concentric shells in each. The 
trajectories of electrons and positrons in the magnetic field were 
simulated directly assuming transversal geometry of the field Bj in 
the jet and B c in the external matter. 

The primary soft photon field (disc and external isotropic) 
is kept constant during the simulations. Any additional soft (syn- 
chrotron) photons produced by the cascade participate in the simu- 
lations in the form of LPs. At the start of simulations, the shell of 
the length Az = 20 between 10 < 2 < 30, and the radial extend 
0.9 < r < 1 is filled by seed isotropic high-energy photons, whose 
energy density is several orders of magnitude less than the energy 
density of the jet. All particles participating in the simulations are 
descenders of these seed photons. In the course of the simulations, 
the jet undergoes differential deceleration. We split the jet into 500 
cylindric shells, calculate the momentum transferred to each shell 
and decelerate each shell independently from others. 
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Figure 6. Total cumulative energy release of the fluid into radiation versus 
time. Solid curve: Example 1; dashed curve: Example 2 (see text). 



5 RESULTS OF SIMULATIONS 

We have made several tens of simulation runs with different pa- 
rameters and various model formulations. In some of them we ob- 
served the exponential energy grow, others gave no effect. As we 
have shown in Section|3| the supercritical behavior appears, when 
the jet Lorentz factor and the density of the soft isotropic radiation 
is sufficiently high. The lowest Lorentz factor, where the cascade 
is possible in principle is around V — 4. Here we present only two 
examples demonstrating the development of the runaway cascade 
for different conditions. 



5.1 Example 1. Weak magnetic field and a "minimal" seed 
radiation 

In this example we assume the disc luminosity Ld,45 = 1, the 
distance of the active region from the black hole Rir = 2, and jet 
parameters V — 10, 9 = 0.05. The total jet kinetic power is equal 
to the disc luminosity r/K — 1 and the Poynting flux is one per 
cent of that 7/b = 0.01. This implies a matter-dominated jet. The 
resulting magnetic field in the jet (comoving frame) is Bj fa 0.35 G 
and we take the external magnetic field B c — 10 -3 G. The ratio of 
isotropic and disc radiation energy densities is ?yi = 0.05, and we 
assume a = 0.4. 

The total (cumulative) energy release as a function of time is 
shown in Fig. [6] We observe a reasonably fast breeding with e- 
folding time t a ~ 0.6. The active layer is rather thin: a half of 
the energy release is concentrated within Sr ~ 0.02 from the jet 
boundary. At t ~ 6 the regime changes: the external shell deceler- 
ates (see Fig.0 and the active layer gets wider (Sr = 0.05 at t = 8 
and Sr — 0.08 at t = 20). The cascade breeding slows down as the 
photon path length through the cycle increases. The total energy re- 
lease into photons reaches 19 per cent of the total jet energy at the 
end of simulation at t = 20. Simulation demonstrated a significant 
pair loading at the late stage: the Thomson optical depth of pro- 
duced pairs across the jet has reached ~ 3 10~ 5 which exceeds the 
initial depth of electrons associated with protons by 25 per cent (as- 
suming proton dominated jet, see Eq. (3|). Most of these pairs are 
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Figure 7. Fluid Lorentz factor versus cylindrical radius of the jet r for Ex- 
ample 1 at various times t = 6 (upper curve), t = 10 (middle curve), 
t = 20 (lower curve). 



loaded in the outer region of the jet (r > 0.8), where they dominate 
the number density of original electrons by a factor of 8. 

While the jet decelerates in our model, the external environ- 
ment is fixed at rest. In reality, it undergoes a radiative acceleration. 
Fig.|8|shows the momentum exchange between matter and photons. 
The momentum transferred to the external environment is an order 
of magnitude less than momentum losses of the jet to radiation. The 
momentum gained by the external shell Sr — 2 10~ 3 around the 
jet boundary is AP ~ 10 48 erg/c. The volume of the innermost 
external shell is 2irRj Sr Az w 10 47 cm 3 . Therefore the volume 
density of the transferred momentum is ~ 10 erg cm _3 /c. If the 
density of external medium exceeds 10 4 cm -3 (i.e. the energy den- 
sity nm p c 2 is higher than 10 erg cm -3 ), we can neglect the effects 
of its acceleration by the deposited momentum. 

The evolution of the photon instant spectrum (i.e. the spec- 
trum of photons which are in the volume at the moment) is shown 
in Fig.|9^. Early spectrum demonstrate two distinct components: 
the TeV Comptonization peak (mainly Compton scattered exter- 
nal isotropic photons) and the synchrotron maximum. After t ~ 6 
the spectrum changes: the main peak moves to lower energies 
and the synchrotron peak declines. The reason for such evolution 
is evident: the system enters a nonlinear stage, because the syn- 
chrotron radiation of the cascade exceeds the initial soft photon 
field. Comptonization losses increase, while the synchrotron losses 
do not change. 

Note, that an observer cannot see the high-energy part of early 
spectra because the isotropic radiation field Fi so is opaque for pho- 
tons with energy e > 1/0. The observed spectrum should have 
such a cutoff in the late spectra. A detection of a sharp cutoff at tens 
of GeV would be the evidence for gamma-rays origin in quasars 
(blazars) at the scale of the broad emission line region. 



5.2 Example 2. Strong magnetic field 

We consider the same parameters as in Example 1, but take r/B 
1 implying magnetically dominated jet with Poynting flux Lb 



Photon breeding in relativistic astrophysical jets 9 



a, 
a, 




0.6 0. 
Radius 

Figure 8. Momentum exchange between the fluid and photons for Example 
1. A narrow distribution (dashes) corresponds to time t = 6, and a wider 
distribution is for t = 20. 
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Figure 9. Evolution of the instantaneous (jet frame) photon spectra within 
the emitting medium, (a) Example 1, (b) Example 2. The observed spectra 
should be close to the latest ones. Time in units Rj /c is marked next to the 
corresponding curves. 



10 45 erg s _1 (for a two-sided jet). The comoving value of magnetic 
field is 3.5 G. The ratio of Compton to synchrotron losses at such 
parameters is low (see Section l3~3l . Therefore, the system remains 
sub-critical and the simulations have shown no supercritical photon 
breeding at such conditions. Our analytical calculations show that 
the maximum of A is about 0.5. The situation changes if we take 
a softer spectrum of the isotropic component with a — 1, then 
max A = 3.5. If we instead increase the soft photon density to 
rji = 0.2, then max A — 1.2. However, the easiest way to achieve 
criticality is to increase the jet Lorentz factor to T = 20, which 



gives maxyl = 8.4. All this options increase the importance of 
Compton cooling relative to the synchrotron (see eg, 1 191 1271 ), 

In the numerical simulations, we followed the last alternative. 
In this case, the comoving magnetic field is 1.75 G. The simulations 
have shown that the active layer is wider than in Example 1: 5r ~ 
0.05 at the beginning of the evolution (therefore this case is more 
stable against the mixing of the boundary layer). The exponential 
growth is slower (see Fig.|6j, but the final energy release is the same 
as in the previous case. The pair loading is an order of magnitude 
smaller than in Example 1 . 

The hard to soft evolution of high-energy peak of the photon 
spectrum shown in Fig.[5J) spans almost all range of peak energies 
observed in blazars. The latest spectrum peaks in MeV range as in 
MeV blazars. 



6 DISCUSSION 

We have demonstrated that a supercritical runaway cascade devel- 
ops under reasonable conditions and can convert at least ~ 20 per 
cent of the jet kinetic energy into radiation. This is certainly not an 
ultimate value: with our simplified model we are able to reproduce 
only the initial stage of the evolution. Indeed, our cylindric shells 
decelerate as a whole with constant Lorentz factor along z axis. 
Therefore, once the outer shells decelerate, the cascade breeding 
slows down everywhere. Actually, the radial gradient of T should 
depend on z and a slow breeding at larger z can coexist with a fast 
breeding at a smaller z. 

The behaviour of the jet is strongly non-linear because the pro- 
cess is very sensitive to a number of details including geometry of 
the magnetic field, density of the external environment, density of 
the isotropic soft photon field, etc. The effect of inhomogeneities, 
which are exponentially amplified, can be dramatic, particularly 
taking a form of flares and moving bright blobs and resulting in for- 
mation of internal shocks. A more realistic model should include a 
detailed treatment of fluid hydrodynamics (at least in 2D) coupled 
with the electromagnetic cascade. 

The model can reproduce the high-energy component of 
blazar radiation. On the other hand, examples presented in this 
work do not reproduce the low-energy synchrotron components as 
prominent as observed in blazars. It is possible that the synchrotron 
bump is produced at much larger distances from the central source 
than the Compton component. 

A certain problem can appear if the jet boundary is turbulent 
(see e.g. lAlov et a l. 1999): the thin active zone at the jet boundary 
layer as in above examples does not exist in this case. Then one still 
can expect to obtain the supercriticality at a moderate transversal 
opacity. However the breeding cycle in this case would take longer 
than Rj/c. If the electromagnetic cascade grows only along the jet, 
the active range of z could be insufficient to provide the growth by 
orders of magnitude. An issue to be studied is whether the cascade 
can grow with time at a fixed z. This could be due to a spatial 
feedback at step 5: a photon from the external environment moves 
upstream to a smaller z than the point where a parent photon has 
been produced at step 3. It is clear that such feedback is weak, but 
in the case of a jet the time is unlimited. Note that even if the su- 
percriticality is not reached, the process can amplify in subcritical 
regime the high-energy output of charged particle acceleration. 

Now let us try to characterize in general terms the mechanism 
we are dealing with. First of all, it belongs to a class of supercriti- 
cal runaway phenomena like neutron breeding in a nuclear pile or 
a nuclear explosion. Such phenomena still seems rather exotic in 
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astrophysics (let alone "trivial" nuclear explosions of supernovae 
or of the accreting matter at the surface of neutron stars). To our 
knowledge there exist only a few works considering such kind of 
phenomena. Stern & Svensson 1 1991) discovered with numerical 
simulations a supercritical behavior of electromagnetic cascade in 
a cloud of ultra-relativistic protons with sufficient compactness. In 
that case the energy is stored in non-radiating protons and the super- 
criticality appears in the energy transfer from protons to pairs and 
photons through photo-meson production. Later th i s mec hanism 
was confirmed analytically by Kirk & Mastichiadis 1 1992). How- 
ever, at that time, there was no clear astrophysical situation provid- 
ing proper conditions. Recently Kazanas et al. (2002) have found 
that a possible site for such phenomenon could be a highly rela- 
tivistic shock in gamma-ra y bursts. 

In this work (see also lSterrl2003l) we propose a different kind 
of a supercritical process where the energy is extracted by particles 
directly from the kinetic energy of the fluid. In principle, the pro- 
to ns can also participate in such mechanism, especially, if a scheme 
of Derishev et al. <2003|) with p + 7 — > n + -k + charge exchange 
works at given conditions. In this case we would have a unified 
mechanism, where the electromagnetic cascade is fed by the fluid 
bulk motion directly and through the high-energy nucleons. 

In any case, supercritical models of energy conversion look 
promising for explanation of such violent phenomena as blazars 
and gamma-ray bursts because the supercriticality does produce vi- 
olent effects. 



7 CONCLUSIONS 

We have proposed and studied a novel photon breeding mechanism 
of the high-energy emission from relativistic jets. We showed that a 
relativistic jet moving through the sufficiently dense soft radiation 
field inevitably undergoes transformation into a luminous state. We 
have considered the application of this mechanism for the AGN 
jets, while actually it may also work in microquasars and gamma- 
ray bursts, if the latter are associated with the well formed narrow 
jets. In general, the mechanism can be characterized as a viscous 
dissipation of the kinetic energy of the jet into high-energy pho- 
tons. We showed that at least 20 per cent of the jet energy can be 
converted into high-energy radiation. From the dynamical point of 
view the mechanism is a supercritical process, which is very similar 
to the chain reaction in the supercritical nuclear pile. The subject to 
exponential breeding in our case is electromagnetic cascade, par- 
ticularly high-energy photons which create a viscosity between the 
jet and the external environment. 

In the case, when the exponential photon breeding does not 
occur, the process still could act as a subcritical amplifier for the 
high-energy output of internal shocks or other mechanisms of par- 
ticle acceleration. Alternatively, it can produce a population of hot 
electrons/positrons which then can participate in the Fermi I type 
of acceleration of charged particles. 

The photon breeding mechanism works very efficiently at the 
following conditions: 

(i) Jet Lorentz factor exceeding 3-4. 

(ii) Weak magnetic field and/or large density of the external soft 
isotropic radiation field. 

(iii) If the magnetic field is strong (equipartition or magnetically 
dominated jet), the supercritical regime requires a dense soft radi- 
ation component (in the jet frame), which requires either a large 
fraction of scattered photons or a large jet Lorentz factor. Alter- 
natively, an additional source of soft photons, e.g. the synchrotron 



radiation from other forms of the jet high-energy activity, e.g. by 
particle acceleration in internal shocks, is needed. The required in- 
tensity of this activity is much smaller than the final high-energy 
output. 

(iv) The presence of a sharp boundary of the jet, with the depth 
of the transition layer less than a few per cent of the jet radius 
(which is still several orders of magnitude larger than required 
for charged particle acceleration). If the jet does not have a sharp 
boundary, then the formation of a supercritical regime is still pos- 
sible at large distances, where the photon-photon pair production 
opacity declines. 
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